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We present the investigation of the neutrino-induced strangeness associated production 
on nuclei in the relativistic plane wave impulse approximation (RPWIA) framework at 
the intermediate neutrino energies. In this study, the elementary hadronic weak am- 
plitudes are embedded inside the nuclear medium for the description of the exclusive 
channels of neutrino-nucleus interactions. These amplitudes are extracted using a model- 
dependent evaluation of the hadronic vertex using the Born term approximation in which 
the application of the Cabibbo V-A theory and SU(3) symmetry are assumed to be valid. 
The nuclear effects are included via the bound state wavefunctions of the nucleon ob- 
tained from the relativistic mean field (RMF) models. Two kinematics settings are used 
to examine various distributions of the differential cross section in the rest frame of 
the target nuclei. The numerical results are obtained for the neutrino-induced charged- 
current (CC) K A-production on bound neutrons in Is^/^ and Ip^^^ orbitals of ^^C. 
The angular distributions are forward peaked under both kinematic settings, whereas 
under the quasifree setting the cross sections tend mimic the missing momentum distri- 
bution of the bound nucleon inside the nucleus. 

1. Introduction 

The study of neutrino interactions, in general, has played a tremendous role in test- 
ing the validity of various models of the electroweak unification theory. In particular, 
the weak neutrino-nucleus interactions involving exclusive associated production of 
strange particles in the intermediated energy region have become one of the major 
theoretical and experimental topics in the fields of astrophysics, cosmology, particle 
and nuclear physics. At the fundamental level of the standard model the theoretical 
and experimental investigations of these exclusive channels assist in gaining new 
and better insight into nucleon and nuclear structures. As such, neutrino interac- 
tions can play a pivotal role by complementing the theories of particle and nuclear 
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physics that have been mainly estabhshed based on the electromagnetic interac- 
tions. Above all, the study of neutrino-induced strangeness associate production 
provides an improved understanding of basic symmetries of the standard model, 
strange-quark content of the nucleon, structure of the weak hadronic form factors, 
strong coupling constants, and medium modification of the elementary amplitudes. 
Therefore, strange particle production via the weak neutrino-nucleus interaction is 
worthy of further theoretical investigation on its own, even more than that due to 
the role it plays as a platform where the studies of particle and nuclear physics can 
overlap. 

Precision neutrino physics necessitates a detailed theoretical knowledge of 
neutrino-induced associated production of strange particles on nuclei owing to the 
presence of such channels as backgrounds in the detection of the neutrino-oscillation 
signal. Along this line MINERj/AIU, K2}M, and MiniBooNEP experiments are re- 
portedly gathering high-statistic data on neutrino-oscillation signals as well as ex- 
clusive kaon-hyperon production cross sections. There are also other relevant back- 
grounds like the quasi-elastic (QE) and pion production channels which have been 
well studied because of their relative simplicity and comparatively bigger cross sec- 
tions. Therefore, these backgrounds need to be carefully addressed in the analysis 
of the neutrino-oscillation data at low and intermediate neutrino energy regime.Sl 
In those experiments the nuclei that have been used as targets are ^^C, ^^O, and 
^°*Pb, just to mention a few. On the other hand, beyond the standard model the 
supersymmetric grand unification theories (SUSY GUTs) clearly emphasizes the 
significance of a thorough understanding of neutrino-induced strangeness associ- 
ated production backgrounds pertaining to the fact that their subsequent decay 
processes, when initiated by atmospheric neutrino fluxes, mimic the nucleon decay 
signal. Currently, there is an ongoing proton decay search being carried out by 
using the world's largest nucleon decay detector Super-KamioKande,^ and future 
proposed plans such as lagunaEI and also other next generation experiments 
involving underground Megaton detectors such as Hyper-KamioKande or UNO.^ 
The success of these experiments depends on how they will be able to separate the 
backgrounds from the decay signal despite the effort to minimize their impact. 

Early experimental observations of the neutrino-induced exclusive final states, 
in which a kaon belonging to a pscudoscalar meson octet is produced in conjunc- 
tion wit h a h yperon in the baryon octet, were reported at AN 
and BNlPl. The data from those experiments, however, are known for their low- 
statistics and large systematic error. Moreover, previous rough theoretical estima- 
tions done by Shroclff^, Mecklenbur JI^, AmeA^, and DewanIS! were based on the 
assumption that nucleon targets are free particles. In actual experiments, however, 
the nuclcons arc bound inside the target, hence invoking the necessity of taking 
into account the medium modification effects such as Fermi motion, Pauli blocking, 
and binding energy correction. As a consequence, a proper theoretical and exper- 
imental co mp arison has never been done owing to those uncertainties incurred on 
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In the intermediate energy region, neutrino-nucleus interactions involves exclu- 
sive channels of either one or two strange particles. The final states of semileptonic 
strangeness production can be categorized in terms of the change in charge and the 
change in strangeness quantum numbers of the weak hadronic transition current. 
The first category comprises the CC and strangeness conserving (AS' = 0) processes 
in which either a K*^ or (S =1) meson is produced in conjunction with either a 
A, E°, or (S = -1) hyperon. The bubble chamber experiments detected events of 
the Vf^n — > /i^K+A channel from this category. The second category is distinguished 
by CC and A5 — 1 currents and here the AS — AQ selection rule is applicable, 
as a consequence only single strange particle or K'^ is produced in the final 
state. Some of the observed channels under this category are z^^p — > /i^K+p and 
Uf^n — > ^~K^p. Neutrino exclusive processes involving the neutral-current (NC) 
and AS* = hadronic weak transition form the third category and the observed 
channels in the bubble chamber experiment are I'^p — > I'^K+A and z/^n — > I'^K'^A. 
So far there has not been any observation of exclusive channels of the forth category 
(i.e., NC and AS* — 1 processes) and henceforth their detection in the ongoing and 
future neutrino experiments will ha ve a great implication for the emergence of new 
physics beyond the standard model.^^^ 

The detailed data with minimum statistical error for heavier nuclei will be avail- 
able from MINERi/A and other ongoing and planned experiments. Moreover, the 
renewed effort of theoretical formulation of the neutrino- induced associated produc- 
tion that incorporates the nuclear effects will allow the comprehensive analysis of 
experimental data usin g rn odels. Motivated by the ever growing experimental inter- 
ests, in our recent worM^ we have revived a few long-forgotten theoretical studies 
of exclusive strangeness associated production via — N interaction with particular 
interest in CC and AS — channels. The paper introduced a new technique of ex- 
tracting weak hadronic elementary amplitudes from a model-dependent evaluation 
of the hadronic vertex in the framework of electroweak theory and fundamental 
symmetries of the standard model. The most general nature of our formulation 
makes it very convenient to be implemented in the studies of the CC and AS" — 1 
and NC and AS = channels as well quite easily, without having to make sig- 
nificant alterations. The model of the current study is constructed on the basis of 
our previous work done at elementary particle level, and aimed at accounting for 
nuclear medium effects. On the other hand, Rafi Alam eta fi^ have investigated the 
CC and AS = 1 processes, in which a single strange particle (5 = 1) is produced, by 
employing the effective Lagrangian formulation in the chiral symmetry framework. 
Therefore, a meaningful interpretation of data and an improved understanding of 
hadronic and nuclear weak physics with particular interest in the Cabibbo V-A 
theory, the SU(3) symmetry, and mcdimn modification effects will be possible in 
the near future. 

It is important to notice the fact that neutrino-nucleus interactions at neutrino 
energy in the order of few GeV also open up to other final state channels, for in- 
stance, a single strange and non-strange meson productions which may even have 
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higher probabihty of detection. Nevertheless, our effort in the theoretical investi- 
gation of the neutrino-induced strangeness associated production of pseudoscalar 
meson and hypcron is motivated and also strongly supported by early experimental 
observations and the present and future experimental activities that we mentioned 
earlier. On the contrary, near threshold, the neutrino-induced associated production 
of vector kaon with identical quark structure as pseudoscalar kaon has never been 
reported, to our knowledge. This may partly be due to their large masses and rel- 
ative instability. Thus the invariant mass of the final state of strange vector meson 
and hyperon requires a minimum neutrino energy that may exceed the near thresh- 
old energy limit upon which our study and those accelerator and non-accelerator 
based experiments are focused. In addition, their instability, leading to very short 
lifetime, makes their experimental observation more difficult than those lighter and 
somewhat stable pseudoscalar kaons. Moreover, at low and intermediate neutrino 
energies nuclear physics studies have greater interest in pseudoscalar mesons than 
vector mesons. From the theoretical perspective, the conventional derivation of the 
invariant matrix element of the associated production has been simplified under the 
psciidoscalar meson field treatment. But for the channels, if any, having the vec- 
tor meson and hyperon fields in the final state, the formulation presents a bigger 
complications and is beyond the scope of this study. 

In this paper a fully relativistic formalism for the production of strange particles 
from nuclei via the weak interaction is presented. The main aim of this formalism is 
to offer a relativistic description in terms of the angular, missing momentum, and 
energy distributions of the differential cross section in the few GeV neutrino energy 
region. We develop a model, that makes use of elementary weak amplitudes, by 
invoking the impulse approximation at the hadronic interaction vertex. We assume 
that the incident neutrino interacts, via the exchange of a single vector boson, 
with a bound nucleon inside the target nucleus and then scatters into a final state 
lepton. On the hadronic part, the electroweak interaction of the weak gauge boson 
with the bound nucleon opens a final state channel involving the production of a 
pseudoscalar meson in association with a hyperon. In this approximation the rest 
of the nucleons behave like spectators and make up the residual nucleus. In this 
specific calculation all external particles will be regarded as free (on-shcll) particles 
and subsequently described by plane wavefunctions except for a bound (off-shell) 
nucleon. The motivation behind this will be expanded upon later. 

In our calculations, the KMF approximation to the Walecka model with the 
NL3 parametrization is employed due to the fact that, unlike a relativistic Fermi 
gas (RFG) model, the nuclear effects such as Fermi motion, Pauli blocking, and 
nuclear binding are naturally incorporated via the nucleon bound state wavefunc- 
tions. On the other hand, the elementary hadronic weak amplitudes are extracted 
from the Born term approximation, in which we assume that the standard Cabibbo 
V-A theory along with the conserved vector current (CVC) hypothesis as well as 
the SU(3) symmetry estimates of the strong coupling constants htc valid and appli- 
cable. In the dynamical description of the weak hadronic vertex, we only consider 
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the lowest-order Feynman diagrams of s-, u- and i-channels mediated by hadrons 
belonging to a baryon octet and a pseudoscalar octet. This approximation is valid 
near threshold where strong resonances are assumed to be absent except for the 
lowest-lying ones, which in fact decay mainly into Ntt and NTTTrffS The numeri- 
cal results of our formalism for a typical neutrino-induced CC reaction process, in 
which K is produced in conjunction with A on a neutron bound in the ^^C nucleus, 
are presented under different kinematic settings. 

This paper is organized as follows. In section [2] we present the relativistic de- 
scription for the neutrino-induced CC associated production on nuclei, whereby we 
incorporated the nuclear structure effects through the bound state wavefunction. In 
section [3] we develop our model for the nuclear weak transition operator in terms of 
the elementary weak current operator by invoking impulse approximation, and then 
evaluate the RPWIA cross section. In section |4] we provide a model-independent 
parametrization of the hadronic current operator in terms of a set of invariant am- 
plitudes and also explain main sources of uncertainty in our formulation as well. In 
order to obtain numerical results we calculate these invariant amplitudes by em- 
ploying the Born term approximation. In section [5] we present the numerical results 
for the RPWIA calculations under different kinematic settings. Finally, in section 
[6] we summarize our results and provide topics for future investigation. 



2. Formalism 

In this section we present the formalism for the relativistic neutrino-nucleus scatting 
processes. The reaction processes under consideration are of the form 

i^ik) + A(P) ^ l{k') + Kip[) + Y{p'^) + B(P'), (1) 

as shown in Fig. [TJ The incoming neutrino with four momentum — {Ek, k) 
scatters into the final lepton I, with k'^ = {Ek', k') off the target nucleus A, with 
= {E, P) via the exchange of a single weak vector boson W = W or Z , that 
carries the four-momentum transfer = (w, q). In the exclusive strangeness asso- 
ciated production channel, the final state involves the production of a pseudoscalar 
meson K, with p'^^ = {Ep[i P'l) accompanied by hyperon Y, with pj'^ = {Ep'^, pj) 
and then both exit the reaction vicinity, leaving behind a recoiling residual nucleus 
B, whose four momentum is labeled as P'^ — {E' , P'). We also denote the masses 
of the target nucleus, final lepton, kaon, hyperon, and residual nucleus as M, , 
M^, My, and M', respectively. 



2.1. Differential cross section 

The theoretical description of the process is given in the laboratory frame in which 
the target nucleus is at rest, and hence its four- momentum becomes P^ = (Af, 0). 
We follow the standard procedure of constructing the differential cross section.i^ 
Thus, for the neutrino-nucleus reaction that proceeds through the KY-production 
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Fig. 1. The lowest order Feynman diagram for the neutrino-induced associated production of 
strange particles on nuclei. The incoming neutrino interacts with the target nucleus via the ex- 
change of a gauge boson, and subsequently the final state lepton, kaon, and hyperon are produced 
and then exit the reaction vicinity, leaving behind the residual nucleus. 



channels as shown in Fig. [T] we may write down the most general differential cross 
section as 

da ^ I ^'^^ I i2ny S^k + P - k' - p[ - p', - P')\M /.p 

|Vl-V2| 

fAr fAr ^ ' 

^(2^)3^(271)3 '(271)3^^(2^)3' 

where |vi — V2I is the relative velocity of the incident neutrino with respect to the 
target nucleus, and in a rest frame of the target nucleus its value becomes unity 
owing to the fact that we neglect the (very small) neutrino mass. The overall four- 
momentum conservation is maintained by the four-dimensional (5-function. A4 ji is 
the invariant matrix element that contains the dynamical information of the weak 
nuclear interaction vertex. Moreover, the expression contains the kinematic phase 
space factor of the final states. Note that Ci are volume normalization factors that 
can be specified according to the Bjorken and Drell convention.!^ 



2.2. Relativistic impulse approximation 

The main focus of this paper will be the derivation of the cross section for the 
associated production of strange particles in the framework of a relativistic impulse 
approximation (RIA). One of the RIA assumptions is that the production process 
takes place via the interaction of a single gauge boson with a bound nucleon inside 
the nucleus as shown in Fig. [21 In addition we make the impulse approximation 
which assumes that the in-medium nuclear current operator can be replaced by the 
free nuclear current operator. It is therefore clear that the model consists of embed- 
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Fig. 2. The vertex approximation of the associated KY-production vertex. Here the vector gauge 
boson interacts with a single bound nucleon inside the target nucleus, and subsequently the final 
state kaon and hyperon are produced, whereas the rest of the nucleons behave as spectators. 

ding the underlying free process in the nuclear medium. However, this immediately 
introduces the problem of nuclear medium effects which must be circumvented by 
making suitable and well-motivated assumptions. The bound state wavefunction of 
the nucleon are obtained from the modified Walecka model by invoking the RMF 
approximation.!^ In general, there are two different frameworks that we can work 
with in the impulse approximation, and in this work we only consider the RPWIA 
that will be discussed in section [31 Note also, that such calculation can be general- 
ized to the relativistic distorted wave impulse approximation (RDWIA) framework, 
which falls beyond the scope of this paper. 

2.3. Kinematics 

It is well established that in the relativistic description of the semileptonic weak 
interaction, the differential cross section can be separated into kinematic and dy- 
namic parts. This, in turn, permits the treatment of these two parts individually 
without loss of generality. In what follows we briefly discuss the degrees of free- 
dom we have at our disposal and the constraints that naturally emerge to minimize 
them. We first define the coordinate system as illustrated in Fig. El The z-axis is set 
along the three-momentum transfer (i.e., z = q). The lepton scattering plane can 
be used to specify the other two coordinate axes: y = k x k / sin^;, where 9i is the 
lepton scattering angle, and consequently the x-axis is constrained to the leptonic 
plane and determined using x = y x z. Recalling that we are working in the rest 
frame of the target nucleus our choice of coordinates gives rise to the vanishing 
j/-components of any three-momenta constrained in the scattering plane. Moreover, 
at the leptonic transition vertex lu ~ Ek — Ek' and q = k — k'. 

The leptonic weak transition presents the challenge of determining the four- 
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Fig. 3. The kinematics of the charge current neutrino-induced associated production in the rest 
frame of the target nucleus. The geometry of the production process is depicted in terms of 
scattering, kaon-production, and hyperon-production planes. For convenience sake the recoil of 
the residual nucleus is left out of this geometric representation of the kinematics. 

momenta fc^, k'^, and given that the kinematic variables E^, w, and |q| are 
fixed. Thus one needs to impose energy-momentum conservation at the leptonic 
vertex and also invoke the on-shell conditions on the neutrino and final state lepton 
in order to determine the three-momenta: 

kx — |k| sin a, ky = 0, kz — |k| cos a (3) 

fc^ = |k'|sin(a-|-6li), fc^ = 0, fc^ = |k'| cos(a + 6*,), (4) 

where 

Ik'l 1 
sin a = -p-j- sin 0; , cosa = — [|k| — |k'| cos6'/], (5) 

and 

Note that we make use of relativistic relations: Ek = |k|, E^, = |k'p + Af^, and 
= — q^. In the Feynman diagrammatic formulation the weak gauge boson 
is identified as an internal line, as a consequence it shall remain off-shell (i.e.. 

The nuclear vertex in Fig. [1] represents the interaction of the weak gauge boson 
with the target nucleus and then resulting in the associated production of strange 
particles. In the final state, in addition to the kaon and hyperon, we also have the 
recoiling residual nucleus. This leads to the general energy-momentum conservation: 

u; + M = Ep,^+Ep,^+E', P' = q-p'i-p^, (7) 

where E' = M' + Tm and we assume that the recoil kinetic energy Tm = p^/2M', 
of the residual nucleus is small enough to be neglected. Note that p,„ = P' is 
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the missing momentum. The mass of the residual nucleus is obtained from M' = 
M — (M„ — El,), where M„ is the nucleon mass and Ei, is the binding energy of 
the nucleon. Note that the binding energy can be determined from different nuclear 
structure models and in this particular work we use the values from the RMF model 
with the NL3 set of parameters. 

Turning our attention to Fig. [2] in which the vertex approximation is applied 
to the weak transition of the nucleus, we notice that the exchange gauge boson 
is absorbed by a single bound nucleon with four- momentum = {Ea, p„), and 
subsequently the kaon and hyperon are produced in their respective planes specified 
by the azimuthal angles (j>[ and 4>2 as shown in Fig.[3l In this reaction picture the rest 
of the nucleons behave as spectators. Furthermore, energy-momentum conservation 
leads to: 



Ea = Ep,^ + Ep>^ - = -P'. (8) 



It is also worth noting that at the nuclear vertex, the target nucleus, residual nu- 
cleus, outgoing kaon and hyperon are on their respective mass shell except for the 
bound nucleon. This is what causes the off-shell ambiguity in connection to ele- 
mentary process. In summary, the choice of a proper reference frame, specifying 
the coordinates, enforcing four-momentum conservation, and invoking the on-shcU 
conditions are used to minimize the number of independent kinematic variables 
needed to describe the reaction process. One also needs to choose a proper orienta- 
tion of the production planes with respect to the scattering plane in order to obtain 
reliable results. 



2.4. Nuclear structure models 

In order to account for the nuclear structure effects, we resort to the RMF ap- 
proximation to the Walecka model. In this model, the bound nucleon is described 
by the exact and self-consistent Dirac-Hartree solutions obtained from the effective 
Lagrangian constructed in the framework of the improved Walecka model. Walecka 
i ntro duced the first quantum hydrodynamics model whi ch i s also known as QHD- 
1.1231 Since then other models such as QHD-lI,l23l25l NL3,'2SI, and FSUGold27 have 
added more complexities to the effective Lagrangian in the RMF approximation to 
improve the model's predictive power, for instance, in reproducing the ground-state 
properties of nuclei with more accuracy. In this work, therefore, we use the NL3 
parameter set to generate the bound state wavefunction of the nucleon, owing to 
the fact that it has been known to be the most successful model in describing the 
ground-state nuclei reasonably well so far. 
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2.4.1. Bound state wavefunctions 



By invoking the RMF approximation to the Walecka model the position space 
bound state wavefunction can be obtained: 



1 



where E is the binding energy, n the generahzed angular momentum, which uniquely 
specifies both the orbital angular monument I, and the total angular momentum j 
where j — \k\ — 1/2; / — n for k > or / = —1 — k, for k < 0. m is the projection 
of J. 3^Km(x) is the spin-angular wavefunction that couples the Pauli spinor Xs n 
with the spherical harmonics Yj^m, (x), of order I: 

3^k™(x) = ^ (/, i,m-s^.,s^.|jm)Y'/,„_s^,(x)Xs^, (10) 

where (Z, ^, m — s^/ , s^/ |j m) are the Clebsch-Gordon Coefficients. Note that g^^ix) 
and fE^(x) are not analytical functions but are generated numerically and then 
must be stored in large arrays. Moreover, the minus sign assigned to k in the lower 
component of the bound state wavefunction has to do with reproducing the exact 
value of the corresponding orbital angular momentum V = 2j — I. 

In the plane wave case (see later) one needs the momentum space bound state 
wave function which is the Fourier transform of the wavefunction in Eq. ((9]) and is 
given by 

Z^..™(P) = / d3xe-'P-ZY,_(x). (11) 

To determine an explicit expression for U^^^ip) we employ the partial wave 
expansion 

e-'P- = 4^ ^ (-1)^ 3, ipx) r,,, (p)r;^ (x), (12) 

L, M 

where (px) is the spherical Bessel function of order L. Therefore, the momentum 
space bound state wavefunction can be obtained: 

UE.mip) = 47r(-z)' , (13) 

V/.JP)'T-P3^_(p)/ 

where the momentum space upper and lower component radial wavefunctions are 
defined as 



9eJp) = J <^^X9EAx)jiipx), (14) 
/e» (P) = sgn(K) / dxx fj^^{x)j^,{px). (15) 
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Figure |4] shows the radial wavefunctions as a function of the magnitude of three- 
momentum for the nucleons in the Is^/^ and Ip^^^ orbitals of ^^C. These wavefunc- 
tions are obtained by using the RMF models mentioned earlier. The plots indicate 
that, under proper kinematic settings, one can explore the missing momentum dis- 
tribution of the cross section and subsequently investigate how the nuclear medium 
effects influence the reaction process. In short, the characteristic curves of the ra- 
dial components wavefunctions suggest that significantly larger cross sections can 
be obtained by carefully tuning the magnitude of the missing momentum in the 
range < 0.3 GeV. 




p (GeV) 



Fig. 4. The momentum space radial wavefunctions of the bound nucleon. The top(bottom) panel 
contains the curves for proton (neutron) bound to the Is^'^^ and Ip^^^ orbitals of ^^C. The 
solid (dashed), dotted(dash-dotted), dash-dot-dotted (short-dashed), and short-dotted (short- 
dash-dotted) curves of gEnifEn) correspond to QHD-I, QHD-II, NL3, and FSUGold parame- 
terizations, respectively, to the RMF model. 
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3. Model for Nuclear Weak Transition 

Apart from RefsP^ ^^^^^I^^^^ there is very little theoretical work on neutrino- 
induced associated production of strange particles. To our knowledge, this paper 
represents the first attempt of a theoretical description of this process from nuclei. 
As such, we develop our model of the nuclear weak transition matrix element for 
KY-production based on the photo- a^^ electroproduction of mesons in association 
with baryons from nuclei .^^^E^EHSIESI in doing so, we invoke the impulse approx- 
imation whereby one assumes that the incident neutrino interacts with a single 
bound (off-shell) nucleon inside nuclear medium via the exchange of a single weak 
gauge boson, and subsequently a pseudoscalar meson is produced in association with 
hyperon. The rest of the nucleons behave like spectators and make up the resid- 
ual nucleus. This assumption leads to the embedding of the elementary hadronic 
weak c urr ent operator^ inside the nuclear medium without changing its on-shell 
details.ISS! Nuclear effects are taken into account using the RMF approximation to 
the Walecka model in the description of the bound state nucleon. 

In similar spirit to the above arguments the conventional interpretation of the 
Feynman diagram representing the vertex approximation in Fig. [5] is necessary. The 
nuclear weak transition matrix element is constructed in terms of the the leptonic 
and nuclear currents, vector gauge boson propagator, and weak vertex coupling 
constants associated to the Feynman diagram that reinforces the assumptions made 
under the impulse approximation. That is, 

-iMf^^i^L^^^H,. (16) 

W 

The leptonic current L^, is completely determined from electroweak theory: 

L^ = u,(fc', /i')7M(^-75)«^(fc), (17) 

where the incident neutrino and the outgoing lepton are described by the helicity 
representation of the Dirac spinor fields. The helicity representation of the Dirac 
spinor of the outgoing lepton can be written as 

This Dirac spinor is non-covariantly normalized and Pauli spinors, 0''~*^(k = 
k(0, (/))), are defined as 

/ cos(6'/2) \ „ /-sin(6'/2)e-*"^\ 

0''-+^(k)= \'[ ], ct>^-\i.)=[ \\\ . (19) 

Vsin(6'/2)e*"*/ V cos(6'/2) / 

Note however, that neutrinos are assumed to be massless, and subsequently they 
can only be described by the negative helicity states. In this paper we do not include 
Coulomb distortions for the CC processes, because of the fact that their effects on 
the reaction cross section are usually assumed to be negligible. 



December 26, 2011 1:58 WSPC/INSTRUCTION FILE ws-ijmpe-adera 



Associated Hyperon- Kaon Production via Neutrino- Nucleus Scattering 13 



In the framework of the impulse approximation, the nuclear weak current is 
written as the sum of the contribution of the elementary hadronic weak c urrents of 
the associated production from the individual nucleons inside the 

nucleusPESl 

H*" = j d3xj^^(x)e*i-^, (20) 

where J'''(x) is the w eak ha dronic transition current for the associated production 
from a single nucleonPSEH 



_(-) 



^'^(x) = (x, Pi) [0'"'(x,p1)J J'^(g)W,_(x), (21) 

where (q) is the elementary hadronic weak current operator that has already been 
developed for the neutrino-induced strange particle production on a free (on-shell) 
nucleonP^ '5'!, \ (x, P2) and cp' ' (x, p[) are the wavefunction of the final state 
hyperon and kaon, respectively, with the outgoing boundary conditions as indicated 
by the minus sign superscripts. At this stage, both of these wavefunctions are kept 
general. Note also that i'2 is an arbitrary spin quantization axis in the rest frame of 
the hyperon and s'2 is the spin projection of the hyperon. The weak coupling factor 
r], may take either of the following values depending on the production channel 
being studied: 

cos0^, for CC& AS = 
% X i sin6'c, for CC&AS = 1 (22) 



Ml ^/2 



i, for NC, 



2 



where = 1.166 x 10 GeV ^ is the Fermi constant, 6^ is the phenomenological 
constant known as the Cabibbo angle (cos^^ « 0.97). Note that this formulation 
is valid in the limit < (Q2 ^ _^2) 



3.1. Relativistic plane wave impulse approximation 

The description of nuclear scattering processes is always challenging due to nu- 
clear effects which must be taken into account in order to establish a reasonable 
theory. Within the framework of the impulse approximation we assume that these 
effects arise when the underlying elementary process is embedded inside the nuclear 
medium. In section [2] we have incorporated the basic effects in terms of the bound 
state wavefunction of the nucleon. There are also other nuclear effects that may 
come into play and their origin can be traced to the final state interactions (FSIs) 
like the strong nuclear force and/or the Coulomb interaction (for CC reactions). In 
this model, however, we assume that the effects form the FSIs are negligible and 
hence the outgoing kaon and hyperon are described by plane waves. 

In this work, we present a formalism based on the RPWIA. This is done for 
a number of reasons. First, it serves as a starting point by providing baseline cal- 
culations against which we can compare the full, if any, RDWIA calculations as 
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well as the high-statistic data that will possibly be available in the near future. 
Second, the inclusion of full distorted waves represents a considerable n umeric al 
challenge as evidenced by similar work done for proton-induced reactions PSESl n 
might therefore be very difficult to disentangle effects related to nuclear structure, 
medium effects, the description of the bound state wavefunction and the elementary 
process as well as distortions all in one go. Finally, despite their significance role 
in the present and future data analysis in the neutrino-oscillation experiments and 
the proton decay searches, the scarcity of theoretical work on the neutrino-induced 
strangeness associated production on nuclei, A{i', IKY)B, is also what has lead us 
to establish this baseline formalism. 

In so doing, we have resorted to using the consequences of a strong analogy 
between the studies of reaction processes that proceed via the weak interactions 
and th at of their counterparts that proceed via the electromagnetic interactions. 
In Refs.^^^^^SEQl^ which the authors extensively investigated the A (7, t:N)B, 
A{e, e'TTN)B, and ^(7, KY)B processes within the RPWIA and RDWIA frame- 
works. What they found out was that the final state distortions only appear to 
suppress the overall magnitude of the RPWIA cross sections without causing any 
substantial change of their shapes. Another example is in Refs.EQSH where a RP- 
WIA formalism was used to investigate the strange quark effects of the axial form- 
factor of the nucleon via neutrino-nucleus scattering. In what follows, therefore, we 
evaluate the invariant matrix element in the RPWIA framework. 

In the plane wave limit the kaon and hyperon are described by the wavefunc- 
tions of the forms: e*Pi '' and (pj, 12, 52)6*^2 respectively, and both are on 
their respective mass shell. Thus in the RPWIA approximation, the nuclear weak 
transition, which is proportional to the Fourier transform of the bound state wave- 
function of the struck nucleon, can be written as 

= U,(p'2, s'^)J^{q)U,_{Yym), (23) 

where p„i = pi + P2 ^ Q and l^EKmiPm) in the momentum space wavefunction 
for bound nucleon as defined in Eq. ()11|) and its explicit form is also given in Eq. 
([T3l) . The weak hadronic current operator, J^(g), is bricfiy presented in section SI 
In the plane wave limit the outgoing hyperon is described by the non-covariantly 
normalized Dirac spinor: 



Uy(p'2,i'2,4) 



where the generalized Pauli spinor with respect to any arbitrary spin quantization 
axis, i = i(a, /?), is defined in terms of the Wigner D-function as 



0^ («,/?)= J2 2?i'(«,/3>0)xJ. 



(25) 
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For the neutrino-induced CC process the outgoing lepton would be a negatively 
charged muon with mass, and hence that makes it detectable by the modern de- 
tectors. The same can be said for the outgoing K -meson. As for the final state 
hyperon, it is only the direction of P2 which is needed to be tracked regardless of its 
magnitude. In contrary, the recoiling residual nucleus goes undetected. Therefore, 
by referring back to Eq. ([2]), we integrate over the recoil three- momentum of the 
residual nucleus, P', and the outgoing hyperon energy, Ep>^, and then the five-fold 
differential cross section can be obtained: 

dV 1 E',E,,E'\k'\\p[\\p',\ 

dEk>dicos9i)dEp,^dnp,dnj,,^ 2(27r)7 |i;' + £;p,(l-p^.(q_p;)/|py2)| I ' 

(26) 

where we have imposed non-covariant normalizations. The norm squared invari- 
ant matrix element is derived from Eq. (|16p . For the KY-production channel that 
proceeds via the CC and AS — process: 

l-^/dcc,...o = ^cos^^,^[L,Lt] ^[^.^.t]^ (27) 

h' m. s'2 

where and i?^ are defined in Eqs. ([T7| and (f23| . In this paper we present results 
for the unpolarized differential cross sections by summing over the initial nucleus 
states, helicity states of final lepton, and spin states of hyperon. 



4. Elementary Hadronic Amplitudes 

The model-independent hadronic weak current operator for the neutrino-induced 
strange particle production from on-shell nucleon is discussed in this section. Note 
that, in this work, we take the same underlying elementary process and embed 
it unmodified inside the nuclear medium (i.e., the impulse approximation). The 
elementary production process of interest is of the form: 

v(k) + N{p) ^ l{k') + K{p',) + Y{p'^). (28) 

The hadronic vertex is represented by the generalized current operator which 
may be constructed based on the global structure the operator processes, i.e., the 
current operator is a Lorentz vector and 4x4 matrix. Having sandwiched the 
current operator between two on-shell Dirac spinors, we can invoke the line ar de- 
composition in the Dirac spinor space in terms of the bilinear covariant 

basisEHl 

/, 7^, 7'^, 7^7'^, (T^^ ■ In doing so, we introduce five unknown vertex amplitudes 
which demand further expansion by using three independent four momenta, the 
metric tensor, and the Levi-Cevita tensor: g^, p^, p'^^j 9^'^ i f^nd e^'^"^ , which we 
have at our disposal. By employing the Dirac algebra and Gordon-like identities on 
the gamma matrices and four vector momenta, the most general representation of 
the weak hadronic current operator takes the following form: 
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+^^75^/+ ^5?'7^ + ^675?'7^, 



(29) 



where 



(30) 



D2q'' + D3P'' + D4P'2'"- 



As similar procedure has already established two such form factors for the two- 
body electromagnetic interactions of Dirac particles. These two independent form 
factors have become a successful means in the study of the structure of nucleons 
in detail by using electron scattering and able to interpret the experimental data 
reasonably well. Our method is also similar in spirit with the CGLN mechani sm 
developed for the photo- and electroproduction of mesons on nucleons and nuclei.S^ 
It is also worth noting that, unlike the electromagnetic production processes, here 
we do not need to impose symmetries such as parity, time-reversal, and charge 
conjugation, since the weak interaction is known to violate each of them. 

4.1. Born term approximation 

The model-independent representation of the weak hadronic current operator for 
KY-production consists of eighteen unknown parametrization amplitudes0 How- 
ever, one needs to resort to model-dependent calculations to determine them. In 
doing so, we employ the Born term approximation, in which the production vertex 
is approximated by a series of the Feynman tree-level diagrams, to describe the 
invariant amplitudes in terms of the hadronic weak form factors and Born cou- 
pling constants. In this approximation we only consider background contributions: 
s-channel via nucleon exchange, u-channel via ground-state hyperon exchange, and 
t-channel via kaon exchange. Note that the labels s, u, and t are the Lorentz in- 
variant kinematic variables: 



The Born diagrams, in general, characterized by three different coupling vertices. 
First, we have a vertex for a single baryon weak transition within the SU(3) baryonic 
octet. Due to this transition three standard form factors: Fi{Q), i^2(Q), and G^{Q) 
come into picture. According to Ref.1321 the vector form factors are determined 
from the corresponding electromagnetic fo rm f actors by using the CVC hypothesis, 
whereas we adopt the method used by Ref.S31 in order to determine the axial form 



s 



{q+pf, t^{q-p[f, 



(31) 



''The detailed derivation of the most general hadronic weak current operator and the extraction 
of these amplitudes from the Born diagrams for CC KY-production processes are given in Ref.H^I 
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factor. The explicit forms of these standard form factors for the case of the weak 
CC transitions amongst the SU(3) baryon octet states are summarized in Table IV 
of RefffS. 

Second, we have the weak transition vertex of a pseudoscalar meson that pri- 
marily comes through the i-channel diagram. This vertex introduces either of the 
following pseudoscalar form factors: Fj^+ (g), -F^+.o (g), or Fj^o,+ (q). We make use of 
the method employed in Rem whereby the author determined their values based 
on their relationship with the vector meson form factors. Finally, we have the strong 
pseudoscalar coupling vertex, that introduces the strong coupling constants, g^vN i 
into the formalism of the process and their values are taken from the broken SU(3) 

symmetry predictions based on the well-known phenomenological coupling constant 
„ 45 46 47 

y-KNN ' 

It is worth mentioning significant uncertainties that our relativistic formalism 
may suffer from. The main uncertainty is presented by the estimated strong cou- 
pling constants used in our calculation owing to the fact that their exact values have 
never been determined to date. The Born approximation is also another source of 
ambiguity due to the lack of improved understanding as to how and which res- 
onance and/or background(non- resonance) terms justifiably contribute to the re- 
action channels of great interest. There are other uncertainties such as non-local 
and off-shell ambiguities incurred by the impulse approximation, specially, as one 
attempts to embed elementary weak hadronic current amplitudes inside nuclear 
medium in which the nucleons are bound. In addition, plane wave limit treatment 
of the final state kaon and hyperon is also a naive approximation in a sense that 
the exclusive cross sections could be sensitive to the distortion effects coming from 
the FSIs. In a future these issues must be addressed in order to give a complete 
theoretical picture of strangeness associated production via the neutrino-nucleus 
weak interaction. 



5. Results and Discussion 

In order to generate numerical results we have chosen a ^^C target and particularly 
investigated the K A-production channel from the bound neutron that proceeds 
via the CC and AS* = process. Note however, that our formalism can easily 
be generalized to other nuclei such as ^He, ^^O, "'"Ca, and ■^°^Pb and production 
channels such as K^E^, K^E , and K E^. The most general form of the differ- 
ential cross section is presented in Eq. ([2]). Nevertheless, one needs to eliminate 
the four-dimensional (5-function which is incorporated to enforce that the overall 
energy-momentum conservation at the nuclear interaction vertex is preserved. Upon 
eliminating the 5-function the five-fold differential cross section is obtained in Eq. 
()26|) . and in what follows this quantity will be used to investigate the neutrino- 
induced associated production process in detail. 

In this paper we present the unpolarized cross section by summing over spin 
states. Note however, that our theoretical framework is also suitable to investigate 
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Fig. 5. Missing momentum distribution of the differential cross section for the quasifree K A- 
production. The short-dashed (short-dash-dotted) line is the contribution from the neutrons in 
Is^/^ {Ip^^'^) orbital of -"^^C and the solid line is their total sum. The kinematics inputs are fixed 
at oj = 1.46 GeV, = 0.05 GeV^, and e{ = 10° . 



polarization observables. As for the spin quantization axis, ij, it is fixed along the 
direction of the three-momentum of the outgoing hyperon (i.e., i'2 = 12(^2; 't'-z))- 
This particular choice is to keep the consistency of our description with regard to 
the weak interaction naturally favoring particles with longitudinal spin polarization 
( i.e., along the direction of its three- momentum), and even the ones with negative 
helicity states in particular to account for the violation of fundamental symmetries 
such as parity by the weak interaction. 

The three-momenta of the outgoing K-meson and A-hyperon are set to lie on 
the same plane ( and for simplicity sake we call it production plane) by imposing 
the constraint: 4>'i = (f>2~ 180 = 4>. Furthermore, we limit ourselves to the coplanar 
setup whereby the scattering and production planes coincide (i.e., (j) — ), since 
this geometrical se tup , as pointed out in the study of photoproduction of associated 
production in Ref.^^, generates larger cross sections relative to the out-of-plane 
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setups. 




Fig. 6. Differential cross section as a function of kaon angle for the quasifree K A-production from 
CC neutrino-nucleus scattering in the intermediate energy region. The kinematics inputs are fixed 
at a; = 1.46 GeV, = 0.05 GeV^, and p,„ = 0.12 GeV. The short-dashed (short-dash-dotted) 
line is the contribution from the neutrons in Is^^'^ (Ip^^"^) orbital of ^^C and the solid line is their 
total contribution. 



It is worthy noting that, the theoretical framework of the associated production 
of strange particles on nuclei presents kinematic flexibility in comparison to the 
underlying production on unb ound n ucleon. In the coplanar setup, we adopt the 
kinematics settings used in First, we implement the setting whereby Ek, 

uj, Q^, O'l, and Pm are specified as input variables. Note that we used Pm in a place 
of |pm|- In this setting one needs to solve for T^^ T^, and For convenience we 
call this arrangement as "setting I" . Second, we make use of the kinematic setting 
that appears to be more natural given the expression of the five-fold differential 
cross section. In this setting Ek, ui, Q"^, T^^, 6[, and 62 are fixed beforehand as 
input variables and then one has to determine Pm and by solving the overall 
four-momentum conservation. This arrangement is labeled as "setting 11" . Under 
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setting II, pm is set to vary in a wider range, and hence we will be able to probe 
the momentum distribution of the bound nucleon. 




Fig. 7. Total differential cross section of quasifree K A-production from CC neutrino-nucleus 
scattering as a function of the missing momentum and outgoing kaon angle at Ej. = 2.5 GeV, 
u) = 1.46 GeV, and = o.05 GeV^. 



5.1. Results of setting I 

This kinematics setting is more like that of the associated production on free nu- 
cleon. Here we investigate the missing momentum and kaon angle distributions of 
the differential cross section for neutrino energies Ek = 1.5, 2.0, 3.0 GeV within 
the plane wave impulse approximation. Under the quasifree treatment we fix w at 
1.46 GeV and at 0.05 GeV^ throughout this section. In Fig. [5] we display the 
missing momentum distribution of the differential cross section of neutrino-induced 
associated K A-production on ^^C at 6[ — 10 . The graphs illustrate the contri- 
butions from neutrons in the Is^^^ and Ip^^^ orbitals of ^^C and their sum for the 

A-production. The results for the "'^^C orbitals clearly justifies the proportional- 
ity between the cross section and the bound state wave function. The contribution 
from the Is^/^ dominates over the Ip^^^ orbital for the missing momentum in the 
range Pm < 0.1 GeV, where as for the rest of the range it is the other way round 
that holds otherwise (up to the overall scaling factor). 

In Fig. [6] we show the quasifree differential cross section as a function of the 
kaon angle and this is done by specifying the magnitude of the missing momentum 
at Pm = 0.12 GeV. The Ip^^^ contribution clearly dominates for the whole range of 
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Fig. 8. Koan angular distribution of the differential cross section for neutrino-induced CC associ- 
ated A-production on ^^C in the intermediate energy region under kinematics setting II . The 
curves correspond to the individual contributions of the orbitals and their sum at oj = 1.32 GeV, 
Q2 = 0.12 GeV^, Tj^ = 300 MeV, and e'^ = 0.5°. 

kaon angle and the curves indicate that the peak occurs at foreword angle. Moreover, 
the curves of the cross section level off slowly suggesting that for the wide range 
of forward angle it is highly likely to detect the outgoing K-meson produced in 
association with A-hyperon on bound neutron struck by incident neutrino with 
intermediate energy. It is important to note that the impulse approximation is also 
more reliable for kaon produced in the direction of forward angles. '^Si4E 

Figure [7] displays the three-dimensional angular and missing momentum distri- 
bution of the differential cross section. This figure provides the wide coverage on 
the characteristic of the differential cross section for the quasifree production of 
strange particles on nuclei. We notice that the larger differential cross section can 
be obtained by tuning the missing momentum in the range Pm < 0.15 GeV and the 
kaon angle in the range 9[ < 30 . The phase space factor suppresses the occurrence 
of the peak at small forward angles and causes the shift in the position of the peak 
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Fig. 9. Angular distribution of the differential cross section with respect to the A-hyperon angle 
for neutrino-induced CC A- production on bound neutron in ^^C nucleus in RPWIA framework 
under kinematic setting 11 at = 5 . 

as opposed to the dynamics factor of the cross section. Nevertheless, both dynamics 
and phase space factors play more or less equal role in determining the shape of 
the cross section for the wider range of the kaon angle. 

5.2. Results of setting II 

Numerical results for the Born term approximation in the plane wave limit are pre- 
sented for the intermediate neutrino energies E/^ — 1.5, 2.0, 2.5, 3.0 GeV. The 
kinematics input values lu — 1.32 GeV and — 0.12 GeV^ are maintained 
throughout the discussion of setting II results. Under this setting we will exten- 
sively explore the angular distributions with respect to kaon and hyperon angels as 
well as energy distributions with respect to kinetic energies of the outgoing kaon 
and muon. 

In Fig. [8] We show the kaon angular distribution for the CC K A-production at 
Tj^ = 300 MeV and ^2 = 0.5 . By implementing the NL3 parametrization scheme to 
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Fig. 10. Differential cross section as a function of kaon and A-hyperon angles at E^. = 2.0 GeV, 
uj = 1.32 GeV, Q2 = o.l2 GeV^, and = 350 MeV. 

a single-particle relativistic shell model, we display the individual contributions of 
the neutrons in the Is^/^ and Ip'^/^ orbitals of "'^^C together with their sum for the 
K A-production. The contributions from the two orbitals are nearly the same up 
to some scaling factor. The distribution curves are peaked in the forward directions 
and they level off faster than the curves for the quasifree production. As pointed out 
in the study of CC neutrino- nucleus scattering,!^ the shapes and peaks of the cross 
sections mirror the missing momentum distribution of the stuck nucleon inside the 
nucleus. 

Similarly, in Fig. |9]we display the A-hyperon angular distribution of the differen- 
tial cross section for neutrino-induced associated production of K A on the bound 
neutron in ^^C target at = 300 MeV and 9[ = 5 . The cross sections from both 
orbitals attain their peak at forward A-hyperon angle indicating the similarity that 
the kaon and hyperon angular distributions may have in common. Moreover, both 
angular distributions fall off sharply before reaching 40 . Furthermore, by tuning 
the kinetic energy and angle of the koan, we are able to shift the peak and even 
obtain a double-peaked A-hyperon angle distributions. 

In order to get a broader insight of our RPWIA calculations under kinematic 
setting II, three-dimensional plots of the differential cross sections are presented 
in Fig. [TOl [121 and [TH These plots represent the total contributions of the ^^C 
orbitals for the CC associated production. Figure [10] displays the three-dimensional 
angular distribution as a function of the kaon and hyperon angles at Ek — 2.0 GeV, 
w = 1.32 GeV, Q2 = o.l2 GeV^, and = 350 MeV. We observe that at fixed value 
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of Tj^ the peculiar shapes and peaks of the cross section are, in fact, in response to 
the missing momentum distribution of the struck neutron in Ip^^^ orbital. 

The richness of kinematics setting II can also be used to shed some light on 
the energy distributions of the differential cross section. Figure [TT] shows the kaon 
kinetic energy distribution of K A-production on bound neutrons in the ^^C orbitals 
at 6[ = 10 and 6*2 = 2.5 . Upon fixing the angles of the outgoing koan and hyperon 
at small forward angles, the kaon energy distribution peaks from the ^^C orbitals 
occur in the range 400 < < 500 MeV. In Fig. [12] we show the differential cross 
section as a function of kaon angle and kinetic energy at 02 = 2.5 . Furthermore, 
Fig. [13] displays the differential cross section as a function of the A- hyperon angle 
and kaon kinetic energy at O'l = 2.5 . These three-dimensional plots provide a 
broader insight for the future experiments to probe the proper kinematics ranges 
that increase the likelihood of detecting the outgoing kaon and hyperon without 
interacting with the residual nucleus. 

Finally, in Fig. [14] we show the muon kinetic energy distribution of the differ- 
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Fig. 12. Differential cross section as a function of kaon angle and kinetic energy. The kinematic 
setting II inputs are fixed at = 2.0 GeV, uj = 1.32 GeV, = 0.12 GeV^, and = 2.5°. 

ential cross section at = 300 MeV and 9[ = 02 = "2.5 . The double-peaked 
structure of the overaU differential cross section can be traced back to the missing 
momentum distribution of the p-shell nucleons of ^^C. We observe that as the muon 
kinetic energy increases, the Is^/^ contribution rapidly levels off and hence leaving 
the Ip"^/^ orbital as the only candidate to extract the interaction information of 
neutrino- nucleus scattering. Note that, a similar remark can also be made for the 
kaon kinetic energy distribution of the differential cross section shown in Fig. 1111 

6. Summary and Conclusion 

We have successfully presented the fully relativistic formalism of the neutrino- 
induced associated production of strange particles on nuclei within the plane wave 
impulse approximation framework. Under this approximation we have assumed that 
the incident neutrino scatters off a single bound nucleon and the remaining nucle- 
ons are treated as spectators. The nuclear medium effects such as Fermi motion, 
Pauli blocking, and nuclear binding are naturally incorporated via the nucleon 
bound state wavefunctions obtained from the improved effective Lagrangian of the 
Walecka model that makes use of the NL3 set of parameters.^^SI 

The leptonic weak current is completely derived in terms of the standard elec- 
troweak theory. However, the same does not apply for the nuclear weak current 
case owing to its complex structure. Therefore, after invoking the impulse approx- 
imation at the interaction vertex, the nuclear current is treated as the sum of 
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the hadronic weak currents, which we re d eveloped for the elementary associated 
strangeness production on free nucleonff^ However, with in our formalism, this ap- 
proximation presents some uncertainties such as off-shell ambiguity, which needs to 
be resolved in a future. In our numerical calculations we have employed the Born 
term approximation at the hadronic interaction vertex. Upon a systematic map- 
ping, the eighteen parametrization amplitudes of the hadronic current operator are 
determined in terms of the standard vector-axial form factors and strong coupling 
constants. 




Fig. 13. Differential cross section as a function of A-hyperon angle and kaon kinetic energy. The 
kinematic setting II inputs are fixed at E^. = 2.0 GeV, lu = 1.32 GoV, = 0.12 GeV^, and 
e' = 2.5°. 



We have shown and discussed our results of the RPWIA calculations under dif- 
ferent kinematic settings in the rest frame of the target nucleus. We investigated 
the angular and energy distributions of the five-fold differential cross section for the 
CC K A-production on bound neutron in ^^C at the intermediate energies. The 
angular distributions peaked at forward angles under both kinematics settings. Un- 
der the first kinematic setting the missing momentum distribution clearly indicates 
the qualitative influence of the bound state wavefunctions have over the Is^/^ and 
Ip"^/^ contributions. On the other hand, the angular distribution shows the dom- 
inance of the p-shell for whole range of kaon angle. This setting does not provide 
a wide kinematic ranges from which the input variables are chosen. Nevertheless, 
we believe that it is a very useful setting in which one can extensively study the 
structure of electroweak interaction and even test the validity of various estimations 
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of strong coupling constants. 
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Fig. 14. Muon energy distribution of the differential cross section for ^^C{u, ^ K A) under 
kinematics setting II in the intermediate energy region. The input variables are fixed at Q2 = 0.12 
GeV2, Tj^ = 300 MeV, and e[ = O'^ = 2.5°. 

The second kinematic setting leaves the missing momentum to vary freely, and 
as a consequence the angular and energy distributions clearly mimic the missing 
momentum distribution of the bound nucleon in the nucleus. This setting offers a 
freedom of tuning the kinematic inputs in the wider ranges and hence causing the 
substantial alteration to the shape and magnitude of the cross sections. Moreover, 
by choosing the input set that favors one orbital contribution over the other, one 
can separately study dynamics of the orbitals. We also believe that this kinematics 
flexibility can offer a wide ground in which one can give a reliable interpretation to 
the experimental data that will be made available in the near future. 

In conclusion, The great merits of the RPWIA formalism are the following: (i) it 
provides a model which incorporates a number of important aspects of the reaction. 
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(ii) it provides a first-order indication of tlie beliavior of tlie cross sections for tfiese 
types of reactions, (iii) from a theoretical perspective it allows a measure of analyt- 
ical manipulation which greatly facilitates the numerical implementation and (iv) 
it will provide the baseline calculations for any possible future model which aims 
to incorporate final-state interactions. Moreover, in the Born approximation of the 
weak hadronic vertex, the inclusion of the relevant resonance contributions in the 
evaluation of the elementary matrix element at the intermediate energy region also 
enhances the importance of our study. One can also take into account the Coulomb 
distortion of the scattered niuon in the CC associated production reactions. Our 
current calculations are timely in a sense that the ongoing neutrino-oscillation ex- 
periments such as MiniBooNE, K2K, and MINERz/A, as well as other experiments 
dedicated for proton decay search like Super-KamioKande, demand an improved 
knowledge of the exclusive reaction channels such as strangeness associated pro- 
ductions. Therefore, our RPWIA results may play a significant role in the analysis 
of data from those experiments, and also in the separation process of backgrounds 
from the proton decay signal. 
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